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ABSTRACT 



The Pierre Auger Observatory has associated a few ultra high energy cosmic 
rays with the direction of Centaurus A. This source has been deeply studied 
in radio, infrared, X-ray and 7-rays (MeV-TeV) because it is the nearest radio- 
loud active galactic nuclei. Its spectral energy distribution or spectrum shows 
two main peaks, the low energy peak, at an energy of 10~^ eV, and the high 
energy peak, at about 150 keV. There is also a faint very high energy (E > 
100 GeV) 7-ray emission fully detected by the High Energy Stereoscopic System 
experiment. In this work we describe the entire spectrum, the two main peaks 
with a Synchrotron/Self-Synchrotron Compton model and, the Very High Energy 
emission with a hadronic model. We consider p7 and pp interactions. For the p7 
interaction, we assume that the target photons are those produced at 150 keV 
in the leptonic processes. On the other hand, for the pp interaction we consider 
as targets the thermal particle densities in the lobes. Requiring a satisfactory 
description of the spectra at very high energies with p7 interaction we obtain 
an excessive luminosity in ultra high energy cosmic rays (even exceeding the 
Eddington luminosity). However, when considering pp interaction to describe 
the 7-spectrum, the obtained number of ultra high energy cosmic rays are in 
agreement with Pierre Auger observations. Moreover, we calculate the possible 
neutrino signal from pp interactions on a Km^ neutrino telescope using Monte 
Carlo simulations. 

Subject headings: Galaxies: active - Galaxies: individual (Centaurus A) - High energy 
cosmic ray: UHECR — radiation mechanism: nonthermal 
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Introduction 



Centaurus A (Cen A) is a Fanaroff & Riley Class I (FRI) active galactic nuclei (AGN). 
At a distance of 3.8 Mpc, it is the nearest radio-loud AGN and an excellent source for 
studying the physics of relativistic outflows and radio lobes. Its giant radio lobes, which 
subtend ~ 10° on the sky, are oriented primarily in the north-south direction. They were 
imaged at 4.8 GHz by the Parkes telescope ( Junkes et al.||1993 ) and studied up to ~ 60 GHz 
by Hardcastle et al. 2009 utilizing the Wilkinson Microwave Anisotropy Probe (WMAP; 



Hinshaw et al. 2009). Cen A has a jet with an axis subtending an angle to the line of sight 



estimated as 15° — 80° (see, e.g. Horiuchi et al. 2006, and reference therein). Cen A has 



been well studied in radio, infrared, optical (Winkler et al. 1975 Mushotzky & Baity 



1976 Bowyer et al. 1970; Baity et al. 1981), X-ray and 7-rays (MeV-TeV) (Hardcastle et 



al. 2003; Sreekumar et al. 1999 Aharonian et al. 2009). A tentative detection (4.5 a) 
of Cen A at very high energy (VHE) in the 1970s was reported by Grindlay et al. 1975. 



Subsequent VHE observations made with Mark HI (Carraminana et al. 1990), JANZOS 



(Allen 1993) and CANGAROO (Rowell et al. 1999; Kabuki et al. 2007) experiments 



resulted in flux upper limits. Cen A was also detected from MeV to GeV energies by all 
instruments on board of the Compton Gamma-Ray Observatory (CGRO) in the period of 
1991-1995, revealing a peak in the spectral energy distribution (SED) in uFi, representation 



at ~ 0.1 MeV with a maximum flux of about ~ 10 erg cm s (Steinle et al. 1998). 



For more than 120 hr Cen A was observed (Aharonian et al. 2005, 2009) by High Energy 



Stereoscopic System (H.E.S.S.) experiment. A signal with a statistical significance of 5.0cr 
was detected from the region including the radio core and the inner kpc jets. The integral 
flux above an energy threshold of ~ 250 GeV was measured to be 0.8% of the Crab Nebula 
(apparent luminosity: L(>250 GeV) ~ 2.6 x 10^^ erg s^^). The spectrum was described by 
a power law with a spectral index of 2.7 ± 0.5stat =t 0.2sys. No significant flux variability 
was detected in the data set. Also, for a period of 10 months, Cen A was monitored by 
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Large Area Telescope (LAT) on board the Fermi Gamma-Ray Space Telescope. Flux levels 
were not significantly different from those found by the Energetic Gamma Ray Experiment 
Telescope (EGRET). However, the LAT spectrum was described with a photon index of 



2.67 ± O.lOstat ± O.OSsys (Adbo et al. 2010). The spectra recorded by the cited gamma-ray 



experiments can be considered as the intrinsic spectra of Gen A because they are not 
affected by the Extragalactic Background Light (EBL) absorption. 

It has been proposed that astrophysical sources accelerating ultra high energy cosmic 
rays (UHEGRs) also could produce high energy 7-rays by proton interactions with photons 
at the source and/or the surrounding radiation and matter. Hence, VHE photons detected 
from Gen A could be the result of hadronic interactions of cosmic rays accelerated by the 



jet with photons radiated inside the jet or protons in the lobes ( Gopal-Krishna et al. 2010 



Rieger & Aharonian 


2009; 


Kachelriess et al. 


2009a|b 


Romero et al. 1996 


Isola et al. 



Honda 


2009 


Adbo et al. 


2010 



Pierre Auger Observatory (PAG) studied the spectra of UHEGR above 57 EeV through 



their shower properties finding a mixed composition of p and Fe (Yamamoto et al. 2007 



Abraham et al. 2008 Unger et al. 2007). By contrast, HiRes data are consistent with a 



dominant proton composition at those energies, but uncertainties in the shower properties 



(Unger et al. 2007) and in the particle physics extrapolated to this extreme energy scale 



(Engel et al. 2007) preclude definitive statements about the composition. At least two 



events of the UHEGRs observed by PAO were detected (Abraham et al. 2007, 2008) inside 
of a 3.1° circle centered at Gen A. 

Synchrotron/Synchrotron-Self Gompton (SSG) models have been very successful in 
explaining the multiwavelength emission from Broad-Line Lacertae (BL Lac) objects 



(Bloom & Marscher 1996 Tavecchio et al. 1998). If FRIs are misaligned BL Lac objects. 



then one would expect synchrotron/SSG to explain their non-thermal spectral energy 
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distribution (SED) as well. In the synchrotron/SSC scenario the low energy emission, 
radio through optical, originates from synchrotron radiation while high energy emission. 
X-rays through VHE 7-rays, originates from SSC However, many blazars have higher 
energy synchrotron peaks, so this mechanism then covers much of the X-ray band; for them 
only the 7-rays come from SSC mechanism. In Cen A, synchrotron/SSC model has been 
applied successfully to fit the two main peaks of the SED, jointly or separated, with one 



or more electron populations (Adbo et al. 2010 Chiaberge et al. 2001 Lenain et al. 



2008 Orellana & Romero 2009; Hardcastle et al. 2009). On the other hand, some authors 



(Dermer et al. 2009 Gupta 2008 Becker & Biermann 2009) have considered hadronic 



processes to explain the VHE photons apparent in the SED. 

In this work we use the fact that leptonic processes are insufficient to explain the 
entire spectrum of Cen A, and introduce hadronic processes that may leave a signature 
in the number of UHECRs observed on Earth. Our contribution is to describe jointly the 
SED of Cen A as well as the observed number of UHECR by PAO. We first require a 
description of the SED up to the highest energies obtaining parameters as: proton spectral 
index proton proportionality constant (Ap) and the normalization energy (Eq). Then, 
we use these parameters to estimate the expected UHECRs observed by PAO. The main 
assumption here, is the continuation of the proton spectrum to ultra high energies. We also 
estimate the neutrino expectation in a hypothetical Km'^ telescope when considering that 
the VHE photons in the SED of Cen A are produced by pp interaction. 



2. UHECRs from Cen A 

The Pierre Auger Observatory, localized in the Mendoza Province of Argentina at 
~ 36° S latitude, determines the arrival directions and energies of UHECRs using four 
fiuorescence telescope arrays and 1600 surface detectors spaced 1.5 km. The large exposure 
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of its ground array, combined with accurate energy and arrival direction measurements, 
calibrated and verified from the hybrid operation with the fluorescence detectors, provides 
an opportunity to explore the spatial correlation between cosmic rays and their sources in 
the sky. The Pierre Auger Collaboration reported an anisotropy in the arrival direction of 



UHECRs (Abraham et al. 2007 2008). While a possible correlation with nearby AGNs 



is still under discussion, it has been pointed out that some of the events can possibly be 



Gorbunov et al. 


2008 


Moskalenko et al. 


2009 


Kachelriess et 



ah 2009a). 



The corrected PAG exposure for a point source is given by 'EtopUj{5s) /^m, where 
Stop = (^) 9 X lO^km^yr, top is the total operational time (from 1st January 2004 until 



August 31st, 2007), uj{5s) — 0.64 is an exposure correction factor for the declination of 



Gen A, and VLqq ~ vr is the Auger acceptance solid angle (Guoco & Hannestad 2008). For 
a proton power law with spectral index Up and proportionality constant Ap, the expected 
number of UHEGRs from Gen A observed by PAG above an energy, Emim is given by. 



^topUj{5s) . f Efnm 

^^UHECR — — -^^TF^ — -C^O 



-Op + l 



(1) 



where Eq is the normalization energy. In other words, the expected number of UHEGRs 
depends on the proton spectrum parameters. If we assume that protons at lower energies 
have hadronic interactions responsible for producing the observed gamma-ray spectra at 
very high energies then we can estimate these parameters. An interesting quantity is the 
apparent isotropic UHEGR luminosity that also depends on the spectrum parameters as. 



Ai^dlApEl f Er, 



(ttp - 2) 



Ea 



2-a„ 



(2) 



where is the distance to Gen A. Gn the other hand, during flaring intervals the apparent 
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isotropic jet power can reach ^ lO^^erg s ^, hence the maximum particle energies of a cold 
relativistic wind with velocity /3, apparent isotropic luminosity (L), Lorentz factor (F) and 



equipartition parameter of the magnetic field (e^) is given by (Dermer et al. 2009) 



Er, 



(3) 



where T = 



3. Leptonic Synchrotron/SSC Model 

In accordance to the AGN model presented in detail by Becker & Biermann 2009, the 
synchrotron photons come from internal shocks in the jet. The associated electrons are 



accelerated by the first Fermi mechanism (Blandford & McKee 1976) and the non-thermal 



electron spectrum can be described by broken power-law given by, 



A. 



dE, " 1 ^-Ee.c 



(4) 



where is the proportionality electron constant, a is the electron spectral index, 
Ee,i = 1e,i^e(? ^ud 7e,i is the elcctrou Lorentz factor. The index ? is m, c or Max for 
minimum, break and Maximum, respectively. For instance, 7e,m is the minimum electron 



^Vietri 


1995; 


Gallant 


1999 



ia - 2) 
1836.15^ -le,T 



[a 



1) 



fes 



548.48(1 + 2) 



1.11 X 10*^(1 + 2 



' obs 



5 X IQ^'^ergs 



-1 



-1/4 



' obs 



dt 



obs 



5 X IQ^^ergs^^ 



2.6 X lO^s. 

1/4 ^ ^^obs 



1/2 



2.6 X lO^s. 



(5) 
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where is the electron energy shock, 5^ = [r(l — is the Doppler factor, = cos^ 

is the observing angle along the line of sight, L°^^ is the observed luminosity, dt°^^ is the 
variability and fes is the ratio of shell expansion time to synchrotron emission time given 



by (Bhattacharjee & Gupta 2003) 



fes{Ee) 



-Be.. 



Ee < E„c 



Ee > Eer 



(6) 



The magnetic field, which comes from an equipartition law, is given by 



B 



L 



obs 



5 X lO'^^ergs' 



-3/2 



dt 



obs 



2.6 X lO^s 



(7) 



As the electrons are accelerated in the shock inside the magnetic field B, they emit 
photons by synchrotron radiation depending on the electron Lorentz factor. The photon 
energy in the source frame is related to the photon energy in the Earth's frame by 



E;'"' = ^E^dDermer & Menon 2009[) then, the observed energies toybicki & Lightman 



1979 ) using equation 5 are given by, 



Eo^s - 12eV^^(l + z)-r^re^e^/^r y/^^it^ 

J^,,m,syn " U- 12 e V _ ^ ^ ^ + <^B^^e^B ^ x 10^3 erg S'W l2.6xl06s. 

f 2 r5 p3 ,-3/2 j obs , _3/2 , Jj-obs 

= 0.01eV(l + z^-5*-^°' ' ^ X / / dt 



■ ' ' /34 V5 X 1043 erg s-V V2.6 x lO^s 

<Ma.,s,n = 4.3xlO«eV(l + z)-i|^ 

We notice that the observed energies correspond to the cut-off energies in the 
synchrotron range and depend on parameters, as fes and 6d, that can be determined 
by fitting the first peak of the SED. On the other hand, the differential spectrum. 



- 9 - 



dN^/dE^, of the synchrotron photons is related to the electron spectrum through 
fes{Ef,) Ee (dNe/dEe) dE^ = E^ (dN^/dE^) dE^, where E^ = C^El and is given in eq. 



(11). Thus, we can obtain the observed synchrotron spectrum as follow (Gupta & Zhang 



2007) 



( rE. 



AM ^obs 

^ dE^J syn 



where 



/ Ej^m,syn \ — 4/3— (a— 3) ( E-y^c.syn N— 1/2/ Ej.syn '>4/3 fpobs ^ rp 

\ En ) \ En ' \ En > ^ ^ ,SVn ^ ^' 



obs 



• E^^c,syn \-l/2( E-i^syn \-{a-3)/2 

' Eo ' ^ Eo ) 

■ E^,syn \-{a-2)l2 

• Eo > 



rpobs ^ Tpobs ^ 

^^,m,syn ^ ^^,syn ^ -^^,c,syn 

'y,c,syn ^ 'y,syn ^ "/,Max,syn 



Ae,^ = 2.16 X 10 



r2 E^ A, e-^--' 



-17 -^0 



1 + ^ 



df 



obs 



d. 



Ce = 1.35 X 10-19 (1 + zf e^i" r-i /_ 

V 5, 



2.6xl06s/ V3.8Mpc 

J^obs 1/2 . fliobs 



(CeEo) 

-1 



a-2 
2 



eV 



-1 



(9) 



(10) 



'11^ 



-5 X 1043ergs-i7 V2.6 x lO^s. 

and is the optical depth (see fig. |3]). Equation [9] represents the low energy contribution 
(IR to optical emission ) to the whole spectrum. 

To obtain the contribution from X-ray to 7-rays, we assume that the relativistic electrons 
inside the jet can upscatter the synchrotron photons in the same knot up to higher energies 
in accordance with, 



2 jpobs 



rpobs 'Z rp< 

'y,m,SSC Je,m ^,m,syni 



2 jpobs 



rpobs 'I 7-11 

-^7,c,S'5C ~ 7e,c ^■y,c,syni 



Jpobs 2 rpobs 

^'y,Max,SSC ~ 7e,Maz ^'y,Max,syn 



:i2) 



Now, with eqs. ([s]) and ([s]) we can finally obtain the inverse Compton photon energies. 



E: 



obs 

'y,m,SSC 



4.01 X 10^ eV ^"" T {1 + z)-' et 6^' 



{a-iy 



L 



obs 



1/2 . (^^obs -1 



5 X lO^Sergs-V V2.6 x lO^s 
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jjiobs 

^■y,c,SSC 



rpobs 

^'y,max,SSC 



M -7/2 rll p7 

3.19 X 10=^ eV (1 + z)-" (1 + x)-^ 



L 



obs 



5.4 X 10^4 eV(l + z)-3^L|Di 



/33 V5 X 10^3 erg s-1 

j^obs X —1/2 ^ dt°'^^ 



-7/2 



dt 



obs 



5 X 10^3 



ergs 



2.6 X lO^s 



2.6 X 106s/ 
(13) 



The Self Synchrotron Compton Spectrum is generally (Fragile et al. 2004) written as, 



A ' ^ 



I — I ~ dEe I \ 

SSC E^ i(j J dEg J \ dEj / syn 



dN^ 



dE^ 



(14) 



Combining Eqs. ^ and ^ with (14), we have that the observed Self-Synchrotron 
Compton spectrum is given by 



E. 



dN^\°bs 



( -S^7, rn. 5SC ^ -4/3- (a-3) ( ^i,c, SSC \-l/2 ( ^i,SSC ^4/3 rpobs ^ rpobs 
y E~o ) y ) I Eo ) ^^,SSC < ^^,m,SSC 



dE^ / SSC 



' Ej^c^SSC \ -l/2( E^, SSC \ -(a-3)/2 
Eo ' y Eo I 

/ -E^7,SSC \-(a-2)/2 
Eo ' 



rpobs ^ rpobs ^ rp 

^-,,m,SSC < -^7,55C < -^7,c,5SC 

rpobs ^ rpobs ^ rpobs 

^T,c,SSC ^ ^T,SSC ^ ^j,rnax,SSC 



(15) 



where. 



A 



■y,SSC 



2.32 X 10"^^ 
,0.505 MeV 



r2 61 Ar^^Al En e-^~' 



/ F N -1/2 

-(q-1) / ^■y,c,syn ^ 



dt 



obs 



dy 



E, 







leV 



2.6xl06s/ V3.8Mpc 



(Ce Eq 



a-2 
I 2 



-'■y,ic 



150 kev) 



(16) 



Summarizing, the Leptonic model describes the whole spectrum at energies below a few 
tens of GeV as stated by equations [9] and 15 considering. 



^]\r -.obs 
^ dE^J Lept model 



^]\J ^obs , AfJ^ohs 



^ / syn 



dE^J SSC 



(17) 
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Hadronic Model 



Some authors (Olinto 2000 Bhattacharjee & Sigl 1998; Stanev 2004 Dermer & 



Menon 2009) have considered possible different mechanisms where protons up to ultra 



high energies can be accelerated. Thus, we suppose that Cen A is capable of accelerating 



protons up to ultra high energies with a power law injection spectrum (Gupta 2008), 



where ap is the proton spectral index and Ap is the proportionality constant. Energetic 



protons in the jet mainly lose energy by p7 and pp interactions (Stecker 1968; Berezinskii 



et al. 1990 Becker & Biermann 2009 Atoyan & Dermer 2003; Dermer & Menon 2009); 



as described in the following subsections. 



4.1. p7 interaction 

The p7; interaction takes place when accelerated protons collide with target photons. 
The single-pion production channels are p + 7 — n + 7r+ and p + 7 — )■ p + vr*^, where 
the relevant pion decay chains are 7r° — )■ 27, tt^ — /x"*" + z/^ — )■ e+ + i^e + ^'/x + ^-'^i and 
TT^ — !■ /i" + z/^ — i- + z/g + t'/x + J^ti (Atoyan & Dermer 2003). 



In this analysis we suppose that protons interact with SSC photons (~ 150 keV) in 



the same knot. If so, the optical depth is given as r. 



p,ssc 



rd6jetTn°l'^^ap^, where is 



the value of the dissipation radius (Bhattacharjee & Gupta 2003), 6 jet is the jet aperture 
angle, = 0.9 mbarn is the cross section for the production of the delta-resonance in 
proton-photon interactions and n^^f^^ is the particle density of SSC photons into the observer 



frame (Becker & Biermann 2009) given by 
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, T obs 

obs _ ('knot 



Assuming that the luminosity of a knot along the jet is a fraction eunot ~ 0.1 of the 
observed luminosity L°^'^ = 5 x 10^^ erg for E'^^ keV, the optical depth is, 



T,,,sc ~ 8.2 X 10-' {^^'^ ( 1^1^°*) ^ io43ergs-0 (lO^^m) (iSOklv) ' 



The energy lost rate due to pion production is (Stecker 1968 Berezinskii et al. 1990), 



^'va = ^ de5^{e)^{e)e dxx~^n{x) (21) 

where n(x) = dn^/de^{e^ = x), <7.„{e) is the cross section of pion production for a 
photon with energy e in the proton rest frame, ^(e) is the average fraction of energy 
transferred to the pion, and eo = 0.15 is the threshold energy, 7p = ep/m^. 

The rate of energy loss, t'^^^, f-nO^p^ ~ ^J'^'p^ (where t'^ ~ r^jV is the expansion time 
scale), can be calculated by following Waxman & Bahcall 1997 formalism. 



obs 



P — P,b 



Here, apeak ~ 5 x 10 cm^ and ^(cpeafe) ~ 0.2 are the values of a and ^ at -E"^ ~ e^eafc 
and Aepeak ~ 0.2 GeV is the peak width. 

The differential spectrum, dN^/dE^ of the photon-pions produced by p7 interaction 
is related to the fraction of energy lost through the equation: f^o(Ep) EpdNp/dEpdEp = 
E^dN^/dE^dE^. If we take into account that vr^ typically carries 20% of the proton's 
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energy and that each produced photon shares the same energy then, we obtain the observed 
gamma spectrum through the following relationship, 



2 dN\ 



TT"— 7 



-1 

Eo J \ Eo 
^ s-ap+2 

I E^ r < E. 



(23) 



where 



Ap^^ = 2.25x10"^^ 



6'^" E^Ap {11. ly-'^^e-^'' 
(1 + zY 



rpobs 



' obs 



150 keV/ V5 X 10^3 ergs- 1 



dt 



obs 



d. 



2.6 X 106s/ V3.8Mpc 
(24) 



and 



212.49 GeV^(^^^ 



(25) 



The eq. 23 could represent the VHE photon contribution to the spectrum. 



4.2. PP interaction 

Hardcastle et al. 2009 argues that the number density of thermal particles within the 
giants lobes is np ~ 10"^ cm"^. If we assume that the accelerated protons collide with 
this thermal particle target then, the energy lost rate due to pion production is given by 



(Atoyan k Dermer 2003), 



where a^p = 30 mbarn is the nuclear interaction cross section, kpp = 0.5 is the inelasticity 
coeficient and n'p is the comoving thermal particle density. The fraction of energy lost by 
PP is f^o^pp ^ t'Jt'pp then. 
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"p ""PP ^ PP 



(27) 



where R is the distance to the lobes from the AGN core. 



The differential spectrum, dN^/dE^ of the photon-pions produced by pp interaction 



is related to the fraction of energy lost through the equation: fnO,pp{Ep) Ep (^^) dEp 



E. 



7 VZeT 



dE^. Taking into account that photon carries 18% of the proton energy, we have 



that the observed pp spectrum is given by (Gupta 2008) 



E' 



dN\ 



r 



dE J 



A 



pp 



ppn 



Ey 



2-ar, 



(28) 



where, 



A. 



pp 



9.97 X 10 



_,,VHlElA,e- 



'l + z) 



R 



100 kpc 



10-4 cm-3 



dt 



obs 



2.6xl06s/ V3.8Mpc 



dy 



(29) 



The eq. 28 could represent the VHE photon contribution to the spectrum. 



5. Calculation of physical parameters and expected UHECRs 



A broadband fit to the SED of Cen A (data from Abdo et al. 2010) using our leptonic 
model (blue line) plus either p7 or pp emission is shown in Figures [T] and [2] respectively. For 
this fit, we have adopted typical values reported in the literature such as luminosity {L°^'^), 



variability {dt ), thermal particle target density (up) and lobes distance (i?)(Adbo et al. 



2010 Dermer et al. 2009 Hardcastle et al. 2009 Romero et al. 1996). The viewing angle 



was chosen in accordance with the observed infrared range (see, e.g. Horiuchi et al. 2006 



and reference therein). Then, from the fit we obtain the values of the bulk Lorentz factor 
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(r), ratio of expansion time (fes), proportionality constants {A^, Aic, Ap^/App), magnetic 
field parameter (^s), electron parameters (^e) and spectral index (a). Other quantities as 
magnetic field {B), electrons Lorentz factors {■je,min, 7e,c ), comoving radius (r^), etc, are 
deduced these parameters. Table 1 shows all the values for used, obtained and deduced 
parameters in and from the fit. 

The fit of the VHE photon spectrum with a hadronic model (either p7 or pp interaction) 
determines the spectral index ctp, energy normalization £'0 and proportionality constant 
Ap (see section 2). So, we calculate the number of UHECRs expected on Earth. Results 
are given in Table 1. As shown, the expected number of UHECRs is extremely high if 
we consider the VHE spectral gamma contribution to come from p7 interactions, while 
considering pp interactions the expected number of UHECRs is in agreement with PAO 
observations. 
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Name 


Symbol 


Value 






Input parameters to the model 


Variability timescales (s) 




2.5 X 10'' 


Adbo et al. 


2010 




Luminosity (erg s"-*-) 




5 X 10"*^ 


Adbo et al. 


2010 





Jet angle (degrees) 6 

Normalization constant (leptonic process) (MeV) Eq 

Normalization constant (p7 process) (TeV) Eq 

Normalization constant (pp process) (TeV) Eq 

Thermal particle target density in lobes (cm"-^) rip 

Lobes distance (kpc) R 



40 
0.1 
1 
1 



Horiuchi et al. 



Jourdain et al. 



Aharonian et al. 



Aharonian et al. 



1 X 10* jHardcastle et al, 

100 



Hardcastle et al. 



Calculated parameters with the model 



Bulk Lorentz factor 


r 


2.06 ± 0.03 


Electron spectral index 


a 


2.837 ± 0.004 


Magnetic field parameter 




0.1073 ± 0.0008 


Electron energy parameter 




0.79 ± 0.14 


Ratio of expansion time 


fes 


0.0385 ± 0.0003 


Proportionality electron constant {eV cm? s)~^ 


A, 


(4.368 ± 0.003) X 10^^ 


Proportionality IC constant (eV cm? s)~^ 




(9.65 ±0.07) X lOif* 


Proportionality proton constant {TeVcm^s)~^ 


A 


(5.9 ±0.4) X 10-'' 


Proportionality proton constant {TeVcm'^s)~^ 




(1.37 ± 0.99) X 10"' 


Proton spectral index 


ap 


2.805 ± 0.008 


Derived quantitatives 


Dopplcr factor 


Sd 


1.47 


Magnetic field (G) 


B 


0.19 


Comoving radius (cm) 


Td 


3.8 X 10^*^ 


Minimum electron Lorentz factor 


7m 


1.36 X 10^ 


Break electron Lorentz factor 


7c 


3.47 X 10^ 


Apparent UHECR Luminosity (p7) (ergs~^) 


Lp 


2.7 X 10"'^ 


Apparent UHECR Luminosity (pp) {ergs^^) 


Lp 


2.9 X 10^^ 


Predicted number of events: p7 interaction 
Predicted number of events: pp interaction 


^ ev,p'y 
^ ev,pp 


8.371 X 10^0 
2.29 



Table 1. Parameters used and obtained from and in the fit of the spectrum of Centaurus A. 
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6. Neutrino expectation for Cen A 

The principal neutrino emission processes in the AGNs are hadronic. These 
interactions produce both, high energy neutrinos and high energy gamma rays, through 
pionic decay. As we mentioned before, hadronic interactions generate mainly pions by 
p + p — 7r° + 7r+ + vr" + X (where X is an hadronic product) and p + 7 — ^ A+ — ?■ tt*^ + 7r+. 
The resulting neutral pion decays into two gamma rays, ttq — )■ 77, and the charged pion into 
leptons and neutrinos, vr^ — )■ + Vfi/v^ + Vfilvy, + z/g/z/e- The effect of neutrino oscillations 



on the expected flux balances the number of neutrinos per flavor (Becker 2008) arriving at 
Earth. Therefore, the measured emission of high energy gamma rays from AGNs suggests 
the possibility to have an equivalent high-energy neutrino flux. In the case of Cen A the 
redshift is 2; = 0.0018, therefore VHE photons are not absorbed from EBL and we can 
consider the observed high energy gamma ray spectra as the intrinsic spectra emitted by 
this source and we can use it for the neutrino flux estimation. 

Concerning the physics environment of Cen A, the optimistic conditions assumed to 
calculate the neutrino expectations are the following, 

1. The high energy gamma ray flux detected by HESS are produced according to the pp 
hadronic scenario in Cen A. 

2. The considered neutrino flux correlated to high energy gamma ray activity has a 
minimum duration of 1 year (i.e. the source is assumed to be stable). 

3. The neutrino spectrum of Cen A is assumed without any cut-off. 

4. The observed gamma-ray spectrum is considered as the intrinsic spectrum of Cen A. 



Considering that neutrinos and gamma rays are produced by the same hadronic 
interaction (pp), we follow the description of Becker (2008) to correlate these two messengers 



- 18 - 



and we assume the neutrino spectrum to be the same as the VHE gamma spectrum recorded 
by H.E.S.S. Therefore we perform a Monte Carlo simulation of a possible Km^ neutrino 
telescope in the Mediterranean sea in order to calculate the expected neutrino event rate. 
We choose this location to have a good sensitivity with respect to the position of Cen A. 

The Monte Carlo simulation takes into account the neutrino source position, the 
propagation of neutrino through water, the charged current interaction with the respective 
muon production, the Cherenkov light produced by the muon, the photons produced by 
the electromagnetic showers and the response of the simulated neutrino telescope. Then we 
calculate the signal to noise ratio in the telescope. 

In this analysis the neutrino "backgrounds" are represented by atmospheric neutrinos 
and cosmic diffuse neutrinos. The atmospheric neutrino "background" is generated by the 
interaction of high energy cosmic rays with nuclei in the atmosphere. The cosmic diffuse 
neutrino "background" is taken as the average rate of neutrinos generated by all the galactic 
and extragalactic non-resolved sources. This cosmic diffuse neutrino flux is discussed by 



Waxman and Bahcall (Bahcall & Waxman 2001 Waxman & Bahcall 1998) and his upper 



limit is given as Eld^/dEy < 4.5 x 10 ^GeVcm ^sr ^. The atmospheric neutrino flux 



implemented in our Monte Carlo is well described by the Bartol model (Barr et al. 2004 



2006) in the range between 10 GeV and 100 TeV. We do not consider the "background" 
from atmospheric muon flux since it is filtered out by the Earth because Cen A is most of 
the time under the horizon for our hypothetical telescope, see FigjSj For the calculation of 
signal to noise ratio we take into account only the "background" inside the portion of the 
sky covered by a cone centered in the Cen A position and having an opening angle of 1° . 
This selection is motivated by the angular resolution of our neutrino telescope. 

Using the assumed neutrino spectrum we obtain for the Km^ telescope the expected 
neutrino event rate shown in Fig. |4j As observed, the integrated signal neutrino event rate 
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in one year of recording data is one order of magnitude below the cosmic neutrino event 
rate and two order of magnitude below the atmospheric neutrino event rate reconstructed 
in the region around Cen A. Moreover, even considering few years of neutrino telescope 
operation, with the considered spectrum, we are not able to disentangle neutrino emission 
from Cen A. 

7. Summciry and conclusions 

We have presented a leptonic and hadronic model to describe the broadband photon 
spectrum of Cen A. Our model has eight free parameters (equipartition magnetic field, 
equipartition electron energy, bulk Lorentz factor, spectral index, ratio of expansion time 
and proportionality constants). The leptonic model describes the spectrum up to a few 
GeV energies while the hadronic model describes the Cen A spectrum at TeV energies. Two 
hadronic interactions have been considered, p7 and pp interactions. In the first case, the 
target is considered as SSC photons with energy of ~ 150 keV, while in the second case, the 
target protons are those in the lobes of Cen A. Only one hadronic interaction is considered 
at the time but in both cases, the proton spectrum is extrapolated up to ultra high energies 
to estimate the number of UHECR events expected at Earth. We have required a good 
description of the photon spectrum to obtain values for the quantities required to estimate 
the UHECR events. 

When p7 interaction is considered, the expected number of UHECR obtained is 
several orders of magnitude above the observed by PAO. However, when pp interaction 
is considered, the expected number of UHECR is in very good agreement with PAO 
observations. 

We have also calculated the neutrino event rate from pp interactions observed by a 
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hypothetical Km^ neutrino telescope in the Mediterranean sea. We have calculated the 
signal to noise ratio considering atmospheric and cosmic neutrino "backgrounds" . We have 
obtained that the expected signal event rate is below the required one to disentangle the 
neutrino emission from Cen A from the " backgrounds" . 

We thank the anonymous referee for the comments given to improve the paper. We 
also thank to Charles Dermer, Markus Botcher, Parisa Roustazadeh, Bin Zhang, Giulia 
DeBonis, Bachir Bouhadef, Mauro Morganti, Dario Grasso, Antonio Stamerra and Teresa 
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This work was supported by DGAPA-UNAM (Mexico) Project Numbers IN112910 and 
IN105211 and Conacyt project number 105033. 



- 21 - 
REFERENCES 

Abraham J. et al. (Pierre Auger Collaboration), 2007, Science, 318, 938 
Abraham J. et al. (Pierre Auger Collaboration), 2008, ApJ, 29, 198. 
Abdo A. A. et al. (FERMI Collaboration), 2010, ApJ, 719, 1433. 
Aharonian F. et al. (HESS Collaboration), 2005, A&A, 441, 465. 
Aharonian F. et al. (HESS Collaboration), 2009, ApJ, 695, L40. 
Allen W.H. et al. (JANZOS Collaboration), 1993, ApJ, 405, 554. 
Atoyan A. M & Dermer C. D., 2003, ApJ, 589, 79. 
Baity W. A. et al., 1981, ApJ244, 429. 

Bahcall J. & Waxman E., 2001, Physical Review D., 64(2): 023002. 

Barr CD., Gaisser T. K. , Lipari P., Robbins S. & Stanev T., Phys. Rev. D 70, 023006 
(2004). 

Barr G.D., Gaisser T. K., Robbins S., & Stanev T., Phys. Rev. D 74, 094009 (2006). 

Becker J. K., 2008, Physics Reports, 458 , 173. 

Becker J. K. & Biermann P. L., 2009, Astropart. Phys. 31, 138. 

Berezinskii, V.S., Bulanov S. V., Dogiel V. A., Ginzburg V. I. & Ptuskin V. S. 1990, 
Astrophysics pf Cosmic Rays (North-Holland: Amsterdam), Ch. 4. 

Bhattacharjee P. & Gupta N., 2003, Astropart. Phys., 20, 169. 

Bhattacharjee P. & Sigl G., 2000, Phys. Rept., 327, 109. 



- 22 - 

Blandford R. D. & McKee C. F., 1976, Phys. Fluids. 19, 1130. 
Bloom, S. D. & Marscher, A. P. 1996, ApJ, 461, 657. 

Bowyer C. S., Lampton M., Mack J. & De Mendonca F., 1970, ApJ161 LI. 

Carramifiana et al., 1990, A&A, 228, 327. 

Cheng K. S. & Wei D. M. 1996, MNRAS, 283, L133. 

Chiaberge M., Capetti A. & Celotti 2001, MNRAS, 324, L33. 

Cuoco A. & Hannestad S., 2008, Phys. Rev. D, 78, 023007. 

Dermer C. D., Razzaque S., Finke J. D. & Atoyan A., 2009, New J. Phys. 11, 065016. 

Dermer C. D. & Menon G. 2009, High energy radiation from Black Holes, Princeton 
University Press. 2009. 

Engel R.l R (The Pierre Auger Collaboration), 2007, arXiv: 0706.1921. 

Fragile P. C, Mathews G., Poirier J. & Totani T., 2004, Astropart. Phys. 20, 591. 

Gallant Y. A., 2002, Lectures Notes in Physics, 589, 24. 

Gopal-Krishna, Biermann P. L. , De Souza V. & Wiitta P. J., 2010, Astropart. Phys. 720, 
L155. 

Gorbunov D., Tinyakov, P. Tkachev, I & Troitsky, S., 2008, Sov. J. Exp. Theor. Phys. Lett, 
87, 461. 

Grindlay J. E, Helmken H. F., Brown R. H.,Davis, J. & Allen L. R., 1975, ApJ, 197, L9. 
Gupta N. & Zhang B. 2007, MNRAS, 380, 78. 
Gupta N., 2008, JCAP, 7060806, 022. 



- 23 - 

Hague J. D., (Pierre Auger Collaboration), 2009, Proc. 31st ICRC, Lodz 

Hardcastle M. J., Worrall D. M., Kraft R. P., Forman W. R., Jones C. & Murray S. S., 
2003, ApJ, 593, 169. 

Hardcastle M. J., Cheung C. C, Feain I. J. & Stawarz L., 2009, MNRAS, 393, 1041. 
Hinshaw et al., 2009, ApJS, 180, 225. 
Honda M., 2009, ApJ, 706, 1517. 

Horiuchi S., Meier D. L., Preston R. A. & Tingay S. J., 2006, PASJ, 58, 211. 
Isola, C, Lemoine M. & Sigl G., 2002, Phys. Rev. D, 65, 023004. 
Junkes E. et al, 1993, ApJ, 412, 586. 

Junkes N., Haynes, R. F., Harnett, J.I., & Jauncy D. L. 1993, A&A, 269,29. 
Kabuki, S. et al. (CANGAROO HI Colaboration) 2007, ApJ, 668, 968. 
Kachelriess M., Ostapchenko S. & Tomas R., 2009, New J. Phys. 11, 065017. 
Kachelriess M., Ostapchenko S. & Tomas R., 2009, Int. J. Mod. Phys. D 18, 1591. 
Lenain J. P., Boisson C, Sol H. & Katarzynski K., 2008 A&A478, 111. 
Moskalenko, I. V., Stawarz, L., Porter T. A. & Cheung, C. C, 2009, ApJ, 693, 1261. 
Mushotzky R. F., Baity W. A., Wheaton W.A., & Peterson L. E., 1976 , ApJ, 206, L45. 
Olinto A. v., 2000, Phys. Rept. 333, 329. 

Orellana M. & Romero G. E., 2009 AIP Conf. Proc, 1123, 242. 

Raue M. & Mazin D. , 2008, International Journal of Modern Physics 17, 1515. 



- 24 - 

Rieger F. M. & Aharonian F. A., 2009, A&A, 506, L41. 

Romero G. E., Combi J. A., Anchordoqui L.A. & Perez S. E., 2009, Astropart. Phys., 5, 
279. 

Rowell, G. P. et al. (CANGAROO III Colaboration) 1999, Astropart. Phys. 11, 217. 

Rybicki G.B. & Lightman A. P. 1979, Radiative processes in Astrophysics (Wiley, New York. 
1979). 

Sreekumar P., Bertsch D. L., Hartman R. C., Nolan P.L. & Thompson D. J., 1999, ApJ, 11, 
221. 

Stanev T., 2004, High Energy Cosmic Rays, Springer, 2004. 
Stecker F. W. 1968, Phys. Rev. Lett., 21, 1016. 
Steinle H. et al., 1998, A&A, 330, 97. 

Tavecchio, F., Maraschi, L. & Ghisellini, G., 1998, ApJ, 509, 608. 

Unger M., Engel R., Schussler F., Ulrich R. & Pierre Auger Collaboration, 2007, Astron. 

Nachrichten, 328, 614. 

Unger M., Dawson B.R., Engel R., Schussler F. & Ulrich R., 2008, Nucl. lustrum. Methods 
Phys. Res. A, 588, 433. 

Vietri M., 1995 ApJ, 453, 883. 

Waxman E. & Bahcall J., 1997, Phys. Rev. Lett. 78, 12. 

Waxman E. & Bahcall J., 1998, Texas simposium on Relativistic Astrophysics and 
Cosmology, December 1998. 

Winkler F. P. & White A. E. 1975 ApJ, 199, L139. 



25 



Yamamoto T. & Pierre Auger Collaboration, 2007, |arXiv: 0707. 2638 



This manuscript was prepared with the AAS lATgX macros v5.2. 



- 26 - 



10"^ E 



111^ I MMII^ 



I MMIII I MMIII I 





I MMIII I MMIII I Minn I 



llJ I 



I I I 



• 


Fermi 


■ 


Swift BAT 




Swift XRT 


• 


Suzal<u 


T 


HESS 




TANAMI 


• 


Arcliival radio-optica 


▲ 


OSSE/COMPTEL 




EGRET 


I 1 


I 1 1 1 1 


10 


10^ 10' 10" 10^ 10 




E(eV) 



I I 1 1 




10* 10"' 10"^ 10"^ 1 



Fig. 1. — Fit of observed spectral energy distribution (SED) of Cen A. The blue line is a fit 
to the broadband SED using the leptonic model described in section 2, while the red curve 
is the p7 contribution described in section 3. 
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Fig. 2. — Fit of observed spectral energy distribution (SED) of Cen A. The blue line is a fit 
to the broadband SED using the leptonic model described in section 2, while the red curve 
is the pp contribution described in section 3. 
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Ey [eV] 

Fig. 3. — Internal optical depth function of photon energy for two different 

luminosities: L°*** = 5.0 x lO^'^erg/s (solid line) and L°^* = 1.0 x lO^^erg/s (dot-dashed line) 
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Km'^S Nu Exp Event Rate [1/year] 




Log^^(E.[GeV]) 

Fig. 4. — Neutrino event rate as a function of neutrino energy. The expected Centaurus A 
event rate for one year of observation with a Km? neutrino telescope assuming the spectrum 
measured by HESS (black line), the event rate for atmospheric neutrinos (green line) and 
cosmic neutrinos (purple line) considering the solid angle of 1° around the sources. A quality 
cut has been applied to reconstructed events. The integrated number of neutrinos is given 
for each case. 
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Fig. 5. — Calculated the Fraction of Day (FoD) when a source with a given declination 
is below the horizon with respect to our hypothetical Km'^ Telescope (solid line) in the 
Mediterranean and ICECUBE experiments (dashed line). For Cen A, FoD is ~ 0.8 when 
Kni'^ Telescope is considered. For ICECUBE, Cen A is all the time above the horizon. 



